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Abstract 

^ ■ This paper addresses the optimal scheduhng of the hquidation of a portfolio using a new 

■ angle. Instead of focusing only on the scheduling aspect like Almgren and Chriss in [2], or 
i only on the liquidity-consuming orders like Obizhaeva and Wang in [33] , we link the optimal 

■ trade-schedule to the price of the limit orders that have to be sent to the limit order book 
to optimally liquidate a portfolio. Most practitioners address these two issues separately: 
they compute an optimal trading curve and they then send orders to the markets to try 
to follow it. The results obtained here solve simultaneously the two problems. As in a 
previous paper that solved the "intra-day market making problem" [2^, the interactions 
of limit orders with the market are modeled via a point process pegged to a diffusive "fair 
price" and a Hamilton-Jacobi-Bellman equation is used to solve the problem involving both 
non-execution risk and price risk. Backtests are carried out to exemplify the use of our 
results, both on long periods of time (for the entire liquidation process) and on slices of 5 
minutes (to follow a given trading curve). 
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in- Introduction 

>: 

Optimal scheduling of large orders in order to control the overall trading costs with a trade- 
^SJ ■ off between market impact (demanding to trade slow) and market risk (urging to trade fast) 

. has been proposed in the litterature in the late nineties mainly by Bertsimas and Lo [8] 

\ and Almgren and Chriss [2]. The original approach has been recently generalized in several 

directions (see for instance [3llll[lll[l5l[l6l[l9l[28l[3Tl[36]); however few attempts have been 
made to drill down the model at the level of the interactions with the order books. The 
more noticeable proposal is the one by Obizhaeva and Wang [33j, followed and generalized 
by Alfonsi, Fruth and Schied [Ij and Predoiu, Shaikhet and Shreve [35]. This branch of the 
optimal trading literatur^ focuses on the dynamics initiated by aggressive orders hitting a 
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\ resilient order-boolU, ignoring trading by passive orders. 
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Recall that during the continuous auction processes implemented by most electronic 
trading pools, market participants send their open interests {i.e. buy or sell orders) to a 
queuing system where a "first in first out" queue stands at each possible price. If a buy 
(respectively sell) order reaches a queue of sell (resp. buy) orders, a transaction occurs (see 
for instance |30j for more explanations and modeling details). Orders generating trades are 
said to be aggressive or liquidity-consuming ones; orders filling queues are said to be passive 
or liquidity-providing ones. In practice, most trading algorithms are as passive as possible 
(a typical balance for a scheduling algorithm is around 60% of passively obtained trades - 
see [29]). 

The economic literature first explored and studied these interactions between orders sent to 
a continuous auction system by different actors from a global efficiency viewpoinlH. 
With the fragmentation of equity markets in the US and in Europe, the issue of linking 
optimal posting prices to the optimal liquidation of a portfolio is more and more important. 
A trading algorithm has to find an optimal scheduling or rhythm for its trading, but also 
to choose a sequel of prices and quantities of orders to send to the markets to follow this 
optimal rhythm as much as possible. 

This paper answers this "optimal scheduling and posting" problem as a whole. Indeed, 
and in contrast with most of the preceding literatur^, we use a new approach, introduced in 
parallel by Bayraktar and Ludkovski [7] in a risk-neutral model, which is liquidity-providing 
oriented: liquidation strategies involve limit orders and not market orders. As a conse- 
quence, the classical trade-off of the literature between market impact, or execution costs, 
and price risk disappears in our setting since no execution cost is incurred. However, since 
the broker does not know when his orders are going to be executed - if at all -, a new risk 
is borne: (non-)execution risk. If a limit ask order is inserted in the order book, probability 
of execution and eventually time of execution will depend on the price of the order. 

In our framework, the fiow of trades "hitting" a passive order at a distance 6 from a 
reference price (the "fair price") St - modeled by a Brownian motion - follows an adapted 
point process of intensity Aexp{—k5). It means that the further away from the "fair price" 
an order is posted, the less transactions it will obtain. In practice, if the limit order price is 
far above the best ask price, the trading gain may be high but execution is far from being 
guaranteed and the broker is exposed to the risk of a price decrease. On the contrary, if the 
limit order price is near the best ask price, or even reduces the market bid-ask spread, gains 
will be small but the probability of execution will be higher, resulting in faster trading and 
less price risk. 

This modeling choice is rooted to a model that has originally been introduced by Ho and 
Stoll [23j and modified by Avellaneda and Stoikov |S] to deal with market making. The main 
characteristic of this model is that it does not explicitly consider the limit order book but 
statistically models liquidity: this is an advantage when it comes to mathematical tractabil- 
ity. This model has been used in the literature on market making by Guilbaud and Pham 
[21], Cartea, Jaimungal and Ricci [IT] and Gueant, Lehalle and Fernandez- Tapia[20j. Our 
paper, along with [3, is an attempt to use it to tackle optimal liquidation issues. However, 
contrary to Bayraktar and Ludkovski [7] who only consider risk- neutral traders, our setting 
takes into account aversion to both price risk and non-execution risk. 

The main limitation of our model is twofold. First, our framework deals with one trading 
venue only. Hence we do not model a "smart order routing" (SOR) mechanism but it can be 
seen as the consolidation of all available trading venues, without taking directly into account 



■^see [T7] for a study of the effect of "smart order routing" on competitive trading venues, or [2] for a study of 
the efficiency of order matching mechanisms. 

^Kratz and Schoneborn [27J proposed an approach with both market orders and access to dark pools. 
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the potential specificity of each of the: 



Second, there is no modehng of the book itself. As a consequence, market orders are ab- 
sent whereas some liquidation strategies should involve both market orders and limit order^. 
Another consequence is that market impact, be it permanent or transient, is not modeled. 
In previous works on optimal trade scheduling, market impact was typically rendered by an 
explicit model because trading was done with liquidity-consuming orders. Similar assump- 
tions can be made on executed orders. However, since the introduction of limit orders in the 
literature on optimal trading is quite recent, there is still no model for the market impact of 
liquidity-providing orders. 

The remainder of this text is organized as follows: in the first part, we present the 
setting of the model and the main hypotheses on execution. The second part is devoted to 
the resolution of the partial differential equations arising from the control problem. Part 
3 deals with three special cases: (i) the time-asymptotic case, (ii) the absence of price risk 
and the risk-neutral case, and (iii) a limiting case in which the trader has a large incentive 
to liquidate before the end. These special cases provide closed-form formulae allowing us 
to better understand the forces at stake. Then, in part 4, we carry out comparative statics 
and discuss the way optimal strategies depend on the model parameters. Finally, in part 5, 
we show how our approach can be used in practice for optimal liquidation, both on a long 
period of time, to solve the entire liquidation problem, and on slices of 5 minutes, when one 
wants to follow a predetermined trading curve. 

1 Setup of the model 

Let us fix a probability space (17, T ^ P) equipped with a filtration {Tt)t>o satisfying the usual 
conditions. We assume that all random variables and stochastic processes are defined on 
(Q,^, (Ji)t>o,P). 

We consider a trader who has to liquidate a portfolio containing a large quantity qq of a 
given stock. We suppose that the reference price of the stock (which can be considered the 
mid-price or the best bid quote for example) moves as a brownian motion with a drift: 

dSt = fidt + adWt 

The trader under consideration will continuously propose an ask quote denoted Sf = St + 
and will hence sell shares according to the rate of arrival of aggressive orders at the prices 
he quotes. 

His inventory q, that is the quantity he holds, is given by qt = Qo — where N"" is the 
jump process counting the number of shares he solqll. We assume that jumps are of unitary 
size and it is important to notice that 1 share may be understood as 1 bunch of shares, each 
bunch being of the same siz^. Arrival rates obviously depend on the price /Sf quoted by 
the trader and we assume that intensity A'* associated to N"" is of the following form: 

A"(5") = Aexpi-kd") = Aexpi-kis" - s)) 

^The actual models studying the SOR problem across several venues are more focused on routing across Dark 
Pools specifically (see [22] for a statistical approach or [51] for a probabilistic one). 

^Limit orders with very low prices are somehow proxies of market orders (see the discussion below). 
^Once the whole portfolio is liquidated, we assume that the trader remains inactive. 

^Typically the average trade size (hereafter ATS) or a fraction of it. If one wants to replace orders of size 1 by 
orders of size Sq in the model, this can be done easily. However, the framework of the model imposes to trade with 
orders of constant size, an hypothesis that is an approximation of reality since orders may in practice be partially 
filled. 
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This means that the cheaper the order price, the faster it will be executed. 

As a consequence of his trades, the trader has an amount of cash whose dynamics is 
given by: 

dXt = {St + 6t)dN^ 

The trader has a time horizon T to liquidate the portfolio and his goal is to optimize 
the expected utility of his P&L at time T. We will focus on CARA utility functions and we 
suppose that the trader optimizes: 



sup E[— exp (—^Xt + qriSr — b))] 

where A is the set of predictable processes on [0,T], bounded from below, where 7 is the 
absolute risk aversion characterizing the trader, where Xt is the amount of cash at time T, 
where qt is the remaining quantity of shares in the inventory at time T and where 6 is a 
cost one has to incur to liquidate the eventual remaining quantitjjfl. 



2 Resolution 

2.1 Hamilton- Jacob i-Bellman equation 

The optimization problem set up in the preceding section can be solved using classical Bell- 
man tools. To this purpose, we introduce the Hamilton-Jacobi-Bellman equation associated 
to the optimization problem, where u is an unknown function that is going to be the value 
function of the control problem: 

(HJB) dtu{t, X, q, s) + ndsu{t, x, q, s) + ^a'^dggu{t, x, q, s) 
+ sup X'^id'') [u{t, x + s + 6'^,q-l,s) - u{t, x, q, s)] = 
with the final condition: 

u{T, x, q,s) = - exp (-7(x + q{s - b))) 
and the boundary condition: 

u{t, X, 0,s) = — exp (— 7x) 

To solve the Hamilton-Jacobi-Belmann equation, we will use a change of variables that 
transforms the PDEs in a system of linear ODEs. 

Proposition 1 (A system of linear ODEs). Let's consider a family of functions {wq)q^fq 
solution of the linear system of ODEs (S) that follows: 

Vg G f^,Wciit) = {aq^ - (3q)wg{t) - r]Wq-i{t) 
with Wq{T) = e~^'^^ and wq = I, where a = ■|7ct^, (3 = kfi and r] = A{1 + ^)~^^'^~\ 
Then u{t, X, q, s) = — ex.p{—^{x + qs))wq{t) is solution of (RJB). 



^We introduced here a constant liquidation cost but most results would be mutatis mutandis the same for an 
instantaneous market impact function b{q). 
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The change of variables used in Proposition 1 is based on two different ideas. First, 
the choice of a CARA utihty function allows to factor out the Mark-to-Market value of the 
portfolio {x + qs). Then, the exponential decay for the intensity allows to introduce Wq{t) 
and to end up with a linear system of ordinary differential equations. 

Now, using this system of ODEs, we can find the optimal quotes through a verification 
theorem: 

Theorem 1 (Verification theorem and optimal quotes). Let's consider the solution w of the 
system (S) of Proposition 1. 

7 

Then, u{t,x,q,s) = — exp(— 7(3; + qs))wq{t) is the value function of the control problem 
and the optimal ask quote can be expressed as: 

2.2 Numerical example 

Proposition 1 and Theorem 1 provide a way to solve the Hamilton-Jacobi-Bellman equation 
and to derive the optimal quotes for a trader willing to liquidate a portfolio. To exemplify 
these results, we compute the optimal quotes when a quantity of shares up to g = 6 has to 
be sold within 5 minutes (Figured]). 




Figure 1: Optimal strategy 6'^*{t,q) (in Ticks) for an agent willing to sell a quantity up to g = 6 within 

5 minutes (^ = (Tick.s~^), a = 0.3 (Tick.s"^), A = 0.1 (s"^), k = 0.3 (Tick"^), 7 = 0.05 (Tick"^) and 

6 = 3 (Tick)) 



We clearly see that the optimal quotes depend on inventory in a monotonic way. In- 
deed, a trader with a lot of shares to liquidate need to trade fast to reduce price risk and 
will therefore propose a low price. On the contrary a trader with only a few shares in his 
portfolio may be willing to benefit from a trading opportunity and will send an order with 
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a higher price because the risk he bears allows him to trade more slowly. 



Now, coming to the time-dependence of the quotes, a trader with a given number of 
shares will, ceteris paribus, lower his quotes as the time horizon gets closer, except near T 
because a certain maximum discount b is guaranteed. At the limit, when t is close to the 
time horizon T, the optimal quotes tend to the same value that depend on the liquidation 
cost 6: (5«*(r,g) = -6+iln(l + |). 

As on the above figure, negative quotes may appear. They appear when the quantity to 
liquidate is large compared to the remaining time, especially when (i) there is a real need to 
liquidate because the liquidation cost b is high and/or (ii) when risk aversion and volatility 
are high, because price risk is then an important consideration. In all cases, it means that 
there is a need to reduce the number of shares at hand and a negative quote should, in 
practice, be seen as an approximation of a market ordei0. 

Also, if we consider the above optimal strategy on a longer time window (see Figure Ej), 
we see that optimal quotes have an asymptotic behavior as the time horizon increases. The 
associated limiting case will be studied in the next section. 




Figure 2: Optimal strategy 6'^*{t,q) (in Ticks) for q = 1,... ,6 and T = 2 hours (/i = (Tick.s 
a = 0.3 (Tick.s-3), A = 0.1 (s^^), k = 0.3 (Tick^^), 7 = 0.05 (Tick^^) and 6 = 3 (Tick)) 

Finally, the average number of shares at each point in time, called trading curve by 
analogy with the deterministic trading curves of Almgren and Chriss, can be obtained by 
Monte-Carlo simulations as exemplified on Figure [3] when the trader uses the optimal strat- 
egy- 



^''in practice, one should not use a liquidation model that uses limit orders when the volume to be executed 
is too large and evidently requires market orders from the very beginning. However, if a negative quote appears 
after some time, then sending a market order to reduce the inventory is a rather acceptable interpretation of the 
model. 
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Figure 3: Trading curve for an agent willing to sell a quantity of shares q = 6 within 5 minutes (/U = 
(Tick.s-i), a = 0.3 (Tick.s-i), A = 0.1 {s'^), k = 0.3 (Tick"^), 7 = 0.05 (Tick"^)) and 6 = 3 (Tick) 



In particular, because 6 = 3, the trader has a weak incentive to liquidate strictly before 
T and there are cases for which liquidation is not complete before time T. This is the reason 
why we do not have E[gr] = on the above figure. We will study below the limiting case 
b —7- +00 that "forces" liquidation before time T. 



3 Special cases 



The above equations can be solved explicitly for w and hence for the optimal quotes using 
the above verification theorem. However, the resulting closed-form expressions are not really 
tractable and do not provide any intuition on the behavior of the optimal quotes. Three 
special cases are now considered for which simpler closed-form formulae can be derived. We 
start with the limiting behavior of the quotes, when the time horizon T tends to infinity. We 
then consider a case in which there is no price risk and a case where the agent is risk-neutral 
to both price risk and non-execution risk. We finally consider, by analogy with the classical 
literature, the behavior of the solution as the liquidation cost b increases. 



3.1 Asymptotic behavior as T — > +00 

We have seen on Figure [2] that the optimal quotes seem to exhibit an asymptotic behavior. 
We are in fact going to prove that 5"'*{0, q) tends to a limit as the time horizon T tends to 
infinity, when the inequality ^ < ^70"^ is satisfieci^. 

Proposition 2 (Asymptotic behavior of the optimal quotes). Let's suppose tha^^fx < ^70"^. 
Let's consider the solution w of the system (5) of Proposition 1. Then: 



hm u;,(0) = ^n — 



r-i>+oo 



9! " - ^ 



^^In particular, when /i = 0, this means that the result is true as soon as ct > 0. 
"'^^This condition is the same as a > /3. 
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The resulting asymptotic behavior for the optimal ask quote of Theorem 1 is: 

lim 5"*(0,g) = iln ^ ^ 



T-5-+00 ' k \ 1 + ? aq^ — /3q 



3.2 Absence of price risk and risk-neutrality 

The above result on asymptotic behavior does not apply to /i = cr = 0. We now concentrate 
on this case in which there is no drift {fi = 0) and no volatility {a = 0). In this case, the 
trader bears no price risk because a = and the only risk he faces is linked to the non- 
execution of his orders. 

We now derive tractable formulae for w and for the optimal quotes: 

Proposition 3 (The no-drift/no- volatility case). Assume that a = and that there is no 
drift (fi = 0). 
Let's define: 

w,{t) = j2t^~""^'~'^iT-ty 

Then w defines a solution of the system (5) and the optimal quote is: 

= _6+iln {l + ^ ^(^-^)' V + 

^ ' k y Y!jll^e-^^^i-i){T -t)i I 7 ^ kJ 

In particular, we see that the asymptotic behavior of Proposition 2 is not true here. We 
also see that the optimal quote is bounded from below by— 6+iln(l + ^). Since execution 
is guaranteed at price s — 6 at time T, it is in particular natural in the absence of price risk, 
that quotes never go below —6. 

Now, if one wants to remove risk aversion with respect to both price risk and non- 
execution risk, one can consider the limit of the above solution when 7 tends to 0. 

One then gets: 

and this is the result of Bayraktar and Ludkovski [7] in the case 6 = 0, because the authors 
do not consider any liquidation cost. 



3.3 Limiting behavior as 6 — +(X) 

Let us now consider the limiting case b — )• +00. Sending b to infinity corresponds to a 
situation in which more and more incentives are given to the trader for complete liquidation 
before time T. If we look at the Almgren-Chriss-like literature on optimal execution, the 
authors are often assuming that qt = Q^i . Hence, if one writes the value functions associated 
to most liquidity-consuming optimal strategies, it turns out that they are equal to —00 at 
the time horizon T except when the inventory is equal to nought (hence b = +00, in our 
framework). However, here, due to the uncertainty on execution, we cannot write a well- 
defined control problem when h is equal to +00. Rather, we are interested in the limiting 



The authors most often consider target problems in which the target can always be attained. 
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behavior when b — )• +00, i.e. when the incentive to Uquidate before time T is large. 



By analogy with the initial literature on optimal liquidation [2], we can also have some 
limiting results on the trading curve. 

Hereafter we denote wi,^g{t) the solution of the system (5) for a given liquidation cost b, 
5^*{t, q) the associated optimal quote and qb^t the resulting process modeling the number of 
stocks in the portfolio. 

Proposition 4 (Form of the solutions, trading intensity and trading curve). The limiting 
solution limfc_>._|_oo Wh^qi^t) is of the form A''Vq{t) where v does not depend on A. 

The limit of the trading intensity \\mi,^^oo Ae^^^b' does not depend on A. 

The trading curve Vb{t) := ^[qb,t] verifies that V{t) = linib^-^-oo Vb{t) is independent of 

A. 

More results can be obtained in the no-volatility case: 

Proposition 5 (no- volatility case, b — )• +00). Assume that a = and consider first the 
no- drift case (fj, = 0). 

hm wb,qit) = "^{T - ty 



\f A \ 



The limit of the optimal quote is given by: 



hm <5-(t,<;)= In — ^-(T-f 
6-i>+oo fc \1 + -j^ q 

The limit of the associated trading curve is V{t) = go (l — ^)^^ 
If H ^ 0, similar results can be obtained: 



^q ( ^p{T-t) 

hm Wb,q(t) = — 

The limit of the optimal quote is: 

hm C ^ = 7 In TV^ p 

The limit of the associated trading curve is V{t) = qo ( ^\^_^-fiT 

4 Comparative statics 

We now discuss the role played by the different parameters on the optimal quotes and illus- 
trate the intuition with actual figures. 



Influence of the drift fi: 
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As far as the drift is concerned, quotes are naturally increasing with ^. If indeed the 
trader expects the price to move down, he is going to send orders at low prices to be executed 
fast and to reduce the impact of the decrease in price on the P&L. This is well exemplified 
by Tabled! 



q 


/i = -0.01 (Tick.s-i) 


fi = (Tick.s-i) 


/i = 0.01 (Tick.s-i) 


1 


9.2252 


10.6095 


12.2329 


2 


6.581 


7.8737 


9.3921 


3 


4.92 


6.1299 


7.5507 


4 


3.6732 


4.8082 


6.1391 


5 


2.6607 


3.728 


4.9765 


6 


1.8012 


2.8073 


3.9806 



Table 1: Dependence on fj, of 6"-*{0,q) with T = 5 (minutes), a = 0.3 (Tick.s 2), A = 0.1 (s k = 
0.3 (Tick-^), 7 = 0.05 (Tick"^) and 6 = 3 (Tick) 

Influence of the volatility a: 



Now, coming to volatility, the optimal quotes depend on cr in a monotonic way. If there is 
an increase in volatility, the price risk increases. In order to reduce this additional price risk 
the trader will send orders at cheaper price. This is what we observe numerically on Table [2j 



q 


cr = (Tick.s 2 ) 


a = 0.3 (Tick.s-^) 


cr = 0.6 (Tick.s 2) 


1 


10.9538 


10.6095 


9.6493 


2 


8.6482 


7.8737 


6.0262 


3 


7.3019 


6.1299 


3.6874 


4 


6.3486 


4.8082 


1.9455 


5 


5.6109 


3.728 


0.55671 


6 


5.0097 


2.8073 


-0.59773 



Table 2: Dependence on a of 6"'*{0,q) with T = 5 (minutes), fi = (Tick.s ^), A = 0.1 (s ^), k = 
0.3 (Tick-i), 7 = 0.05 (Tick-^) and 6 = 3 (Tick) 



Influence of the intensity shape parameter k: 



Now, as far as k is concerned, the dependence of the optimal quote on k is ambiguous 
because the interpretation of k depends on the optimal quote itself. An increase in k cor- 
responds indeed to a decrease in the probability to be executed in most cases the model is 
used, that is for quoted prices above the reference price. However, due to the exponential 
form for the execution intensity, the exact opposite is true for quoted prices below reference 
price, an increase in k implying an increase in the probability to be executed for quotes 
below reference price (see the above discussion on negative 6 in section 2.2). 

In the case of positive optimal quotes, an increase in k then forces the trader to decrease 
the value of the optimal quotes as exemplified by Table [3l However, if price risk is really 
important (high volatility and/or large quantity to liquidate) the optimal quotes may be 
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negative and, in that case, the role of k is reversed. This is the case when a takes (unreal- 
istically) high values, as on Tabled] for q >2. 



q 


k = 0.2 (Tick ') 


k = 0.3 (Tick 


k = 0.4 (Tick ') 


1 


15.8107 


10.6095 


7.941 


2 


11.9076 


7.8737 


5.7972 


3 


9.4656 


6.1299 


4.4144 


4 


7.6334 


4.8082 


3.3618 


5 


6.1436 


3.728 


2.5011 


6 


4.8761 


2.8073 


1.7688 



Table 3: Dependence on k of 6^*{0,q) with T = 5 (minutes), fi = (Tick.s ^), a = 0.3 (Tick. 
A = 0.1 (s-^), 7 = 0.05 (Tick-i) and 6 = 3 (Tick) 



q 


k = 0.2 (Tick"') 


k = 0.3 (Tick"') 


k = 0.4 (Tick"') 


1 


2.8768 


0.79631 


-0.031056 


2 


-4.0547 


-3.8247 


-3.4968 


3 


-8.1093 


-6.5278 


-5.5241 


4 


-10.9861 


-8.4457 


-6.9625 


5 


-13.2176 


-9.9333 


-8.0782 


6 


-15.0408 


-11.1488 


-8.9899 



Table 4: Dependence on k of 6"'*{0,q) with T = 5 (minutes), ^ = (Tick.s ^), cr = 3 (Tick. 
A = 0.1 (s-^), 7 = 0.05 (Tick"i) and 6 = 3 (Tick) 

Influence of the intensity scale parameter A: 



Now, coming to A, we observe numerically, that the optimal quote is an increasing func- 
tion of A (see Table E]). If A increases, the probability to be executed indeed increases and 
the trader will then increase his quotes to obtain transactions at higher prices. 



q 


A = 0.05 (s-i) 


A = 0.1 (s-i) 


A = 0.15 (s-i) 


1 


8.4128 


10.6095 


11.9222 


2 


5.6704 


7.8737 


9.1898 


3 


3.9199 


6.1299 


7.4491 


4 


2.5917 


4.8082 


6.1302 


5 


1.5051 


3.728 


5.0525 


6 


0.57851 


2.8073 


4.1341 



Table 5: Dependence on A of (5"*(0,g) with T = 5 (minutes), = (Tick.s"^), a = 0.3 (Tick. 
k = 0.3 (Tick"^), 7 = 0.05 (Tick~^) and 6 = 3 (Tick) 

Influence of the risk aversion j: 



Turning to the risk aversion parameter 7, two effects are at stake that go in the same 
direction. The risk aversion is indeed common for both price risk and non-execution risk. 
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Hence if risk aversion increases, the trader will try to reduce both price risk and non- 
execution risk, thus selling at cheaper price. We indeed see on Table [6] that optimal quotes 
are decreasing in 7. 



q 


7 = 0.01 (Tick"^) 


7 = 0.05 (Tick~^) 


7 = 0.5 (Tick-^) 


1 


11.2809 


10.6095 


9.84 


2 


8.8826 


7.8737 


6.7461 


3 


7.4447 


6.1299 


4.7262 


4 


6.4008 


4.8082 


3.189 


5 


5.5735 


3.728 


1.9384 


6 


4.8835 


2.8073 


0.88139 



Table 6: Dependence on 7 of 6''*{0,q) with T = 5 (minutes), /i = (Tick.s-^), a = 3 (Tick.s'a), 
A = 0.1 (s-^), k = 0.3 (Tick"^), and 6 = 3 (Tick) 

Influence of the liquidation cost b: 



Finally, the influence of the liquidation cost b is straightforward. If b increases, then the 
need to sell strictly before time T is increased because the value of any remaining share at 
time T decreases. Hence, the optimal quotes must be decreasing in b and this is what we 
observe on Table \7\ 



q 


6 = (Tick) 


b = 3 (Tick) 


6 = 20 (Tick) 


1 


10.7743 


10.6095 


10.4924 


2 


8.0304 


7.8737 


7.7685 


3 


6.278 


6.1299 


6.0353 


4 


4.9477 


4.8082 


4.7229 


5 


3.859 


3.728 


3.6509 


6 


2.9301 


2.8073 


2.7374 



Table 7: Dependence on b of 5'^*{0,q) with T = 5 (minutes), /U = (Tick.s ^), cr = 3 (Tick.s 2)^ 
A = 0.1 (s-^), k = 0.3 (Tick-^) and 7 = 0.05 (Tick^^) 



5 Applications 

Before using the above model in reality, we need to discuss some features of the model that 
need to be adapted before any backtest is possible. 

First of all, the model is continuous in both time and space while the real control prob- 
lem under scrutiny is intrinsically discrete in space, because of the tick size, and in time, 
because orders have a certain priority and changing position too often reduces the actual 
chance to be reached by a market order. Hence, the model has to be reinterpreted in a 
discrete way. In terms of prices, quotes must not be between two ticks and we decided to 
round the optimal quotes to the nearest ticlJ^. In terms of time, an order is sent to the 
market and is not canceled nor modified for a given period of time At, unless a trade occurs 



^^We also, alternatively, randomized the choice with probabilities that depend on the respective proximity to 
the neighboring quotes. 
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and, though perhaps partially, fills the order. Now, when a trade occurs and changes the 
inventory or when an order stayed in the order book for longer than At, then the optimal 
quote is updated and, if necessary, a new order is inserted. 

Now, concerning the parameters, a, A and k can be calibrated on trade-by-trade limit 
order book data while 7 has to be chosen. However, it is well known by practitioners that 
A and k have to depend at least on the actual market bid-ask spread. Since we do not 
explicitly take into account the underlying market, there is no market bid-ask spread in the 
model. Thus, we simply chose to calibratJ^ A and k as functions of the market bid-ask 
spread, making then an off-model hypothesis. 

As far as 7 is concerned, a choice based on a Value at Risk limit is possible but requires the 
use of Monte-Carlo simulations. We decided in our backtests to assign 7 a value that makes 
the first quote 5""* equal to 1 for average values of A and k. 

Turning to the backtests, they were carried out with trade-by-trade data and we assumed 
that our orders were entirely filled when a trade occurred at or above the ask price quoted by 
the agent. Our goal here is just to provide examples in various situations and, to exemplify 
the practical use of this model, we carried out several backtestJ^ on the French stock AXA, 
either on very short periods (slices of 5 minutes) or on slightly longer periods of a few hours. 
We believe that the model is particularly suited to optimize liquidation within slices of a 
global trading curve, be it a TWAP, a VWAP, or an Implementation Shortfall trading curve. 

The first two examples (Figures H] and [5]) consist in liquidating a quantity of shares equal 
to 3 times the AT^^. The periods have been chosen to capture the behavior in both bullish 
and bearish markets. 

On Figure HI we see that the first order is executed after 50 seconds. Then, since the 
trader has only 2 times the ATS left in his inventory, he sends an order at a higher price. 
Since the market price moves up, the second order is executed in the next 30 seconds, in 
advance on the average schedule. This is the reason why the trader places a new order 
far above the best ask. Since this order is not executed within the time window At, it is 
canceled and new orders are successively inserted with lower prices. The last trade happens 
less than 1 minute before the end of the period. Overall, on this example, the strategy works 
far better than a market order (even ignoring execution costs). 

On Figure O we see the use of the strategy in a bearish period. The first order is executed 
rapidly and since the market price goes down, the trader's last orders are only executed at 
the end of the period when prices of orders are lowered substantially since selling becomes of 
utmost importance. Practically, this obviously raises the question of linking a trend detector 
to these optimal liquidation algorithms. 



^^A and k are updated as new information comes. 

^^No drift in prices is assumed in the strategy used for backtcsting. 

^^In the backtests we do not deal with quantity and priority issues in the order books and supposed that our 
orders were always entirely filled. 
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14:00 
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Figure 4: Backtest example on AXA (November 5*'* 2010). The strategy is used to sell a quantity of 
shares equal to 3 times the ATS within 5 minutes. Top: quotes of the trader (bold line), market best 
bid and ask quotes (thin lines). Trades are represented by dots. Bottom left: evolution of the inventory. 
Bottom right: cash at hand. 
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13.76 



13.75 



13.74 - 



13.73 



13.72 



13.71 



13.7 



13.69 - 



13.68 




15:00 



15:01 



15:02 



15:03 



15:04 



15:05 



Figure 5: Backtest example on AXA (November 5*'* 2010). The strategy is used to sell a quantity of 
shares equal to 3 times the ATS within 5 minutes. Top: quotes of the trader (bold line), market best 
bid and ask quotes (thin lines). Trades are represented by dots. Bottom left: evolution of the inventory. 
Bottom right: cash at hand. 
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Finally, the model can also be used on longer periods and we exhibit the use of the 
algorithm on a period of two hours, to sell a quantity of shares equal to 20 times the ATS 
(Figure [6]) . 



13.7 - 
13.68 - 




14:00 14:30 15:00 15:30 16:00 14:00 14:30 15:00 15:30 16:00 



Figure 6: Backtest example on AXA (November 8 2010). The strategy is used to sell a quantity of 
shares equal to 20 times the ATS within 2 hours. Top: quotes of the trader (bold line), market best bid 
and ask quotes (thin lines). Trades are represented by dots. Bottom left: evolution of the inventory. 
Bottom right: cash at hand. 
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Conclusion 



As claimed in the introduction, this paper is, to authors' knowledge, the first proposal to 
optimize the trade scheduling of large orders with small passive orders when price risk and 
non-execution risk are taken into account. The classical approach to optimal liquidation, 
following the Almgren-Chriss framework, consisted in a trade-off between price risk and 
execution cost/market impact. In the case of liquidity-providing orders, this trade-off dis- 
appears but a new risk is borne by the agent: non-execution risk. 

The problem is then a new stochastic control problem and an innovative change of vari- 
ables allows to reduce the 4-variable Hamilton-Jacobi-Belmman equation to a system of 
linear ordinary differential equations. Practically, the optimal quote can therefore be found 
in two steps: (1) solve a linear system of ODEs, (2) deduce the optimal price of the order 
to be sent to the market. 

We studied various limiting cases that allowed to find the asymptotic behavior of the 
optimal strategy and to find the result obtained in parallel by Bayraktar and Ludkovski [7], 
taking the risk- neutral limit. This also allowed us to confirm our intuition about the role 
played by the parameters. 

Numerical experiments and backtests have been carried out and the results are promis- 
ing. However, two possible improvements are worth the discussion. 

First, no explicit model of what could be called "passive market impact" {i.e. the per- 
turbations of the price formation process by liquidity provision) is used here. Just note that 
up to now, no quantitative model of this type of impact has been proposed in the litera- 
ture. Thanks to very promising and recent studies of the multi-dimensional point processes 
governing the arrival of orders (see for instance the link between the imbalance in the order 
flow and the moves of the price studied in [T2] or [TB], or interesting properties of Hawkes- 
like models in [6]), we can hope for obtaining models of this kind in the near future. The 
authors will try to embed them into the HJB framework used here. In between, an on-going 
work dedicated to model dependencies between the Brownian motion supporting St and the 
Poisson process A'^" is under consideration. 

Second, the separation of the variables (x, s) and (t, q) is a property associated to the use 
of a CARA utility function and to the brownian dynamic of the price and is independent of 
the exponential decay for the arrival of orders. An on-going work aims at generalizing the 
above model to a general function A"(-) using this separation of variables. 
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Appendix 



Proof of Proposition 1 and Theorem 1: 

First, let us remark that a solution {wq)q of (S) exists and is unique and that, by im- 
mediate induction, its components are strictly positive for all times. Then, let's introduce 
u{t, X, q,s) = — exp {—j{x + qs)) Wq(t)~ k . 



We have: 



dtu + jidgU + -a dggU 



l Wq{t) -i^(J^ 2 

- — j-u - jq/iu + — — g u 

kWq{t) 2 



Now, concerning the non-local part of the equation, we have: 

sup X"-{6^) [u{t, X + s + 6"',q - l,s) - u{t, x, q, s)] 

ga 



= supyle '^^'^ u{t,x,q^s) 



The first order condition of this problem corresponds to a maximum and writes: 

(, + ,)exp(-,n(^^)"*=* 
Hence we introduce the candidate 5"* for the optimal control: 



Wq{t) 



+ 



and 



k \Wq-i{t)J 7 



supA''((5") [u{t,x + s + 6"-,q - l,s) - u{t,x,q,s)] 



7 



= ~A 



fc + 7 

7 



A exp{-kd"-*)u{t, X, q, s) 



Hence, putting the three terms together we get: 



dtu{t, X, q, s) + ndsu{t, x, q, s) + -a dggu{t, x, q, s) 



+ sup X"-{6"-) [u{t, X + s + 5", g - 1, s) - u{t, x, q, s)] 

ga 



7^(7(*) 7^cr^ 2 A 7 



-J Wq-i{t) 

Wg{t) 



u 



7 u 



k Wq(t) 



Wq{t) + kliqWq{t) ^—q^Wq{t)+A[l + ^j ^ ^^Wg-l(t) 







Now, noticing that the boundary and terminal conditions for Wq are consistent with the 
conditions on u, we get that u verifies (HJB). 



Now, we need to verify that u is indeed the value function associated to the problem 
and to prove that our candidate is indeed the optimal control. To that purpose, let's 

consider a control u ^ A and let's consider the following processes for r G [t,T]: 



18 



dS'/ = ^idT + adWr, Si'' = s 



dX. 



t,x,u 



{Sr + Vr)dN^, X, 



t'XM 



where the point process has stochastic intensity (At-)t with A,- = Ae~^'^^l^ ^il^. 
Now, let us write Ito's formula for u since u is smooth: 

u{T, Xtf^_f ,qrj?^ , S!^') = u{t, X, q, s) 

+ jj"^ (n(T, X*'!'^ + S'f + z.., q'ff - 1, S'/) - u{t, S'/)) Kdr 

+ j\dsu{T,X';l'\q';'h\S'/)dWr 

+ (n(r, X*'!'^ + 5*'^ + q'ff - 1, - u{t, X'f_:\q'f_l\ 5*'^)) dM" 
where M"" is the compensated process associated to N"" for the intensity process {\t)t 



Now, we have to ensure that the last two integrals consist of martingales so that their 
mean is 0. To that purpose, let us notice that dsU = —^qu and hence we just have to prove 
that: 



E 



■ pT 

J uir,X^'J ^q;"!! ,5"; ) dr 



< +00 



E 



and 



E 



I I vt,x,i/ t,q,u c<t,s\ 

u{r,X;_' ,q;'!! ,5"; ) 



X-j-dr 



Xrdr 



< +00 



< +00 



Now, remember that the process q^'i''^ takes values between and q and that t G [0,T]. 
Hence, 3e > 0, Wq{t) > e for the values of t and q under scrutiny and: 



< e exp (^-27(x - q\W |U + 2<7^mf^^ 'li„f^,(,,^, 5*-<o; 

< exp {-2j{x - g||z^-||oo)) (l + exp (^-2^q ^mf^^ 5*'^^^ 



Hence: 



E 



Jt J Ut 



<e-^exp(-27(x-g||z.-||oo))(T-t) ( 1 + E 



exp ( -27g inf 5*'" 



This intensity being bounded since v is bounded from below. 
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<e-^exp(-27(x-g||z^-||oo))(r-t) ( 1 + E 



exp —2^q inf S ' 
Te[t,T] 



< e"^ exp (-27(3; - g||z/-||oo))(T - t) (l + e~2'^'?'E exp (27go-VT-t|y 



< +00 

where the last inequahties come from the reflection principle with Y ~ AA(0, 1) and the fact 
that E [e'^l^l] < +00 for any C e M. 

Now, the same argument works for the second and third integrals, noticing that v is 
bounded from below and that A is bounded. 



Hence, since we have, by constructio: 



2 

a / x^t.a;,!/ t,q,v ot,s\ , o / vt,x,u t,q,u cit,s\ , a2 i vt.x.u t.q,u c<t,s\ 



u{t, Xl'^'" + Sl'^ + z.*, g;-!''^ - 1, S'/) - n(r, q'^^ , 5*'^) ) A. < 



we obtain that 



E 



frp -Yt^x^u t,q,i/ Qt^s^ 



E 



u[l,X^J ,qrfz ) '-'t j <u[t,x,q,s 



and this is true for all G Since for = (5"* we have an equality in the above inequality 
we obtain that: 



sup E 



ii(T X^'"^'" n^''^'" < ii(f r n - F \v(T Y^''''^'"' n*'^'^"* S"*'"*^ 

This proves that u is the value function and that 6'^* is optimal. 
Proof of Proposition 2: 
We have that 



□ 



\fq e fi,Wg{t) = {aq^ - l3q)wg{t) - r]Wg-i{t) 



Hence if we consider for a given Q S N the vector w{t) 

w'{t) = Mw{t) where: 

/ 

-T] a-P ■•. 



(wo{t)\ 
wi{t) 



we have that 



M 



\ 







/ 1 \ 



with 'w{T) 



-kb 



-T] a(g-l)2-/3(Q-l) 
-r] aQ'^-PQj 



. Hence we know that, if we consider a basis (/o, • • • , /q) of eigen- 



vectors {fj being associated to the eigenvalue aj^ — there exists (cq, . . . ,cq) G M'^+i 
independent of T such that: 



"'^^This inequality is also true when the portfolio is empty because of the boundary conditions. 
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j=0 

Consequently, since we assumed that a > P, we have that w°° := Um7-_^+oo w{0) = cofg. 
Now, is characterized by: 

{aq^ -Pq)w^=VW^_„q>0 < = 1 

As a consequence we have: 

The resulting asymptotic behavior for the optimal ask quote is: 



lim ,5«*(0,g) = ^In ' ^ ^ 



T-^+oo ^ ' ^ k \1 + ^ aq"^ - /3q ^ 

□ 

Proof of Proposition 3: 

The result of Proposition 3 is obtained by induction. For q = the result is obvious. 
Now, if the result is true for some q we have that: 



Hence: 



9+1 , 



J- 



9+1 , 



j=0 



This proves the results for w and then the result follows for the optimal quote. □ 
Proof of Proposition 4: 

Because the solutions depend continuously on h, we can directly get interested in the 
limiting equation: 

Mq e N, Wq{t) = {aq^ - I3q)wq{t) - ■qwq-i{t) 
with Wq{T) = lg=o and w;o = 1- 

Then, if we define Vq{t) = limb_^_)_oo -^jr^, v solves: 

Vg G n,Vq{t) = {aq'' - Pq)Vq{t) - fjvg^i{t) 
with Vq{T) = lq=o and vq = 1, where 77 = ^ is independent of A. 
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Hence Vq{t) is independent of A. 
Now, for the trading intensity we have: 



lim Aexpi-k6i*{t,q)) = lim ^^^^^^ f i + 



vM) V ^ 



7 

(t) V ' k 



and this does not depend on A. □ 
Proof of Proposition 5: 



Using the preceding proposition, we can now reason in terms of v and look for a solution 



of the form Vq{t) — — 



Then, 

"iq en,Vq{t) = -PqVg{t) - fiVq-i{t), Vq{T) = lq=o, Vq = 1 

^ h'{t) = -ph{t) - fi h{T) = 

Hence, if /3 = fc^ / 0, the solution writes Vq{t) = |!.(£iPM^z^))i:i)9. 
From Theorem 1, we obtain the limit of the optimal quote: 

b^+oo \k \q p J J \ k 

Using the expression for ?), this can also be written: 

1 / A l e^(^-*)-l \ 



k 

Now, the for the trading intensity we get: 



lim Aexp(-Mr(t,g)) = (I + t)?- ^ 



Hence, because the limit of the intensity is proportional to q, the limit V{t) of the trading 
curve is characterized by the following ODE: 

Solving this equation, we get: 

V{t) = goexp(^-(l + ^)^ \,,^£, J s 



pef^ 1 

/ /, . 7\ , /, J- \ 

goexp 



goexp 



1 

?yJe'3(T-t)^ 



\ 14-^ 

l_e-^(^-*)\ 



When /? = (ie. /X = 0) we proceed in the same way or by a continuity argument. □ 
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